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ABSTRACT 


We have tested the implications and limitations of Program 
ACRETE, a scheme based on Newtonian Physics and accretion with 
unit sticking efficiency, devised by Dole (1970) to simulate the 
origin of the planets. The dependence of the results on a variety 
of radial and vertical density distribution laws, on the ratio of 
gas to dust in the solar nebula, on the total nebular mass, and on 
the orbital eccentricity, e, of the accreting grains are explored. 

Only for a small subset of conceivable cases are planetary systems 
closely like our own generated. Many models have tendencies towards 
one of two preferred configurations: multiple star systems, or 
planetary systems in which Jovian planets either have substantially 
smaller masses than in our system or are absent altogether. But for 
a wide range of cases recognizable planetary systems are generated — 
ranging from multiple star systems with accompanying planets, to 
systems with Jovian planets at several hundred AU, to single stars 
surrounded only by asteroids. No terrestrial planets were generated 
more massive than 5 Earth masses. The number of planets per system 
is for most cases of order 10, and, roughly, inversely proportional 
to e. All systems generated obey a relation of the Tltius-Bode 
variety for relative planetary spacing. The ease with which planetary 
systems are generated, using such elementary physical assumptions, 
supports the idea of abundant and morphologically diverse planetary 
systems throughout the Galaxy. 



I . INTRODUCTION 


The ultimate problem In planetary studies Is the origin 
of the solar system. Despite a serious recent attack on the problem 
by many Investigators — much of which has been published in 
the pages of Icarus over the last few years — it seems safe 
to say that no generally acceptable detailed model of the origin 
of the solar system exists. Indeed, the rate of change of. models 
of origins, even in the hands of experienced individual investigators 
is a clear indication of 'the uncertainty of the 
subject. Furthermore, almost all of the detailed models have 
concentrated on the important early stages of solar system history, 
particularly the solar nebula* and not on the origins of planets 
per se . For example, the significant paper by Goldreich and Ward 
(1973) carries the history of the solar nebula up to the generation 
of planetesimals of about the size of Phobos and Deimos. On the 
observational side, the most recent work has served to cast doubts 
on the reliability of claimed identifications of extrasolar planetary 
systems (Gatewood, 1976), rather than providing a data source 
against which models of origins can be tested. Under these circumstances 
any model which purports to generate planetary systems recognizably 
similar to our own deserves careful scrutiny. 

Any acceptable model for the formation of the solar system 
should be able to account at least for its most obvious characteristics: 
the distinction between terrestrial and Jovian planets, the spacing 
of planetary orbits and the distribution of planetary mass with 
heliocentric distance. Such parameters as the rotation periods of the 
planets; the orbital configuration of comets, asteroids and the particles 
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In the rings of Saturn; and the anomalous obliquities of Uranus and Venus 
are presumably details not essential to an understanding of the 
formation processes — although it is possible that they might 
provide significant clues. 

To the best of our knowledge, the only existing model 
wh'Ch attempts to generate mature planetary systems as opposed 
to other stages in the evolution of solar nebulae is that of 
Dole (1970). In his computer simulation, "accretion nuclei” 
of specified mass eire injected in prograde orbits in the invariable 
plane of a primitive solar nebula composed of both gas and dust. 

The physics are simply Newtonian mechanics and perfectly inelastic 
collisions. When accretion nuclei collide with dust grains the 
grains adhere with unit efficiency. Growing accretion nuclei 
beyond a certain mass gravitationally accrete gas as well. When 
two accretion nuclei collide they stick also and produce a larger 
planetesimal. The process is permitted to continue until all the 
dust and some of the gas is gathered into planets. For some choices of 
input parameters the resulting planetary configurations (Figure 1 ) 
are remarkably like those of the solar system. 

At first sight it appears quite extraordinary that so simple 
a physical protocol can lead to so recognizable a set of planetary 
systems. The program takes no explicit account of chemical segregation 
with heliocentric distance in the solar nebula, of a clearing out 
by solar radiation pressure and the solar wind of the inner solar 
system during the T Tauri stage of the sun, of hydromagnetic effects, 
turbulent convection, or of planets 

dynamically unstable because of rapid rotation. If the origin of 
the planets can indeed be understood with such elementary assumptions 
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and with plausible Input parameters, the model deserves much 
deeper attention. 

Dole was able to produce planetary systems of recognizable 
characteristics only with a certain choice of Input parameters 
and assumed structure of the solar nebula. The present paper is 
devoted to a critical examinations of these assumptions and an 
exploration of the consequences of variations of parameters and 
assumed solar nebular structure. 

For example, the coplanar character of the simulated 
planetary systems is a direct consequence of the fact that the 
accretion nuclei are Injected with zero Inclination in prograde 
orbits. The formation and dynamical properties of the accretion 
nuclei are not further justified by Dole; we will 
discuss them further in the light of more recent research. The 
increased current skepticism (Gatewood, 1976) ''on earlier 
reductions of perturbations in the proper motion of Barnard's 
Star removes the props from the reduction by Black and Suffolk 
(1973) according to which the planets of the Barnard Star system 
would not have been in coplanar orbits. 

II. THE COMPUTER MODEL 

The computer simulation program, called ACRETE, was written 
by J. Rice and generously provided to us by S. Dole. We, have varied 
the program where necessary. In this section we describe the 
essential features of ACRETE. 

l.The solar nebula is taken to have the shape of an "exocone" ^ 
seen edge-on in Fig. 2. The shape is assumed to arise from an 
originally spherical cloud of gas and dust with some nonzero net 
angular momentum in which dust particles with bits highly inclined 



to the Invariable plane are eventually degraded to orbits of 
lower inclination through inelastic collisions. Most models of 
the solar nebula assume either a similar configuration or a 
cylindrical (disk-shaped) distribution of matter in which the 
density of gas and dust falls off away from the central plane. 

In its original application, ACRETE does not take account of the 
vertical density distribution. A corrected treatment will be discussed 
in a later section. 

2. The mass ratio of gas to dust in the nebula Is a constant, 

K = pg / Pd where pg and p^j are respectively the radially-dependent 
density of gas and dust. While it may be reasonable to expect that 
this ratio will be Independent of radial distance in the central 
plane of the nebula, the mass difference between a dust particle and 
a gas molecule will insure different scale heights for their respective 
vertical density distributions. This correction, however, would 
complicate the computer program greatly, and was not taken, into 
account either by Dole or by us. Dole used the value K = 50, which 
will be shown later to be a reasonable number. 

3. In Dole’s study, the density distribution of dust is 

»/ * 

p^ = Ae where r is the distance from the center of mass of 

the cloud in astronomical units (AU), and A and a are adjustable 
parameters. In addition to experimenting with A and a, we have 
varied the functional form itself. Dole was able to generate 
aesthetically pleasing (i.e. solar system-like) planetary systems 
when A = 0.0015 Mq/AU^ and a = 5. The justification for these 
particular choices was one of convenience; Dole was not striving for 
any generality, since he stated that the object of his exercise 
was to generate planetary systems similar to our ov;n. We will 
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examine other choices of A and o as well as other choices for 
the functional form p^. 

The dust particles comprising the cloud (other than 
the accretion nuclei) are all given the same orbital eccentricity 
e, an input parameter (Dole's value is e = 0,25)» and are 
taken to have randomly distributed semimajor axes and inclinations. 

We will examine the consequences of other choices of e, but for 
simplicity will not assume a distribution function for various 
values of the orbital eccentricity of dust grains. 

5. The accretion nuclei are taken to have some initial mass 

m which is an input pai-ameter of the program. The nuclei are 
0 

injected into prograde orbits of zero Inclination, with semimajor 
axes randomly distributed between 0.3 and 50 AU, and with eccentricities 
given by the distribution function e = l-(l-Y)*^*®^^, where Y is 
random between zero and one. This form is an empirical distribution 
derived by Dole which reproduces the distribution of planetary 
eccentricities in the solar system. The small exponent yields small 
eccentricities and, since any nucleus undergoing accretion would 
suffer numerous inelastic collisions that would tend to circularize 
its orbit, it was not felt necessary to change the exponent even though 
it is an input parameter. Similarly, the bounds of 0.3 and 50 AU for 
the semimajor axes are also input variables, but changing them changes 
neither the physics of the problem, nor (to any substantial degree) 
the results. Occasionally, these limits were moved closer together 
for convenience when no planets could be formed at the extremities 
of the cloud. 

A nucleus captures all dust particles which cross its orbit 
(sticking coefficient unity), plus those whose orbits fall in 
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an unstable region related to its gravitational cross section 

The radial extent of this region around the accretion nucleus 

1 /•* 

is given by x “ ru > where r is the distance of the nucleus 

from the center ot the nebula and y is its reduced mass with 

respect to the Sun: u = m/(l + m), where 'm is the nuclear mass 

expressed in solar masses. The expression for x is an approximation 

to the,-.solution of the restricted three-body problem. Birn (19’’3) 

finds the exponent to be 1/3 instead of 1/4, but this was not 

changed in the program, since the effect of the change can be 

shown to be small. Also, it is implicitly assumed that the semimp.Jor 

axes of all orbits precess through all directions in the 

Invariable plane via accumulated gravitational perturbations. 

7. Nuclei accrete only dust initially, until their masses (and 

hence escape velocities ) are high enough to permit the retention 

of gas as well. If we assume that an accreting planetoid of mass m 

has uniform density, its escape velocity Vg is proportional to m^/3. 

A gas molecule at temperature T has a velocity “ and, if we 

— 1 /P 

assume a temperature-distance dependence of T(r) = (r/r^)” ' 

(where r is the radial distance from the central scar), then the 
gas velocity becomes v The functional form chosen for T (r) 

O 

is appropriate for an optically thin solar nebula and for some choices 
of optically thick nebulae. For retention of gas above some critical 
mass m^j, we demand v^ >v .» or » so that m^, = , 

where C]_, € 2 % are constants of proportionality. practice, 
is 10“5 when m^. is measured in solar masses and 

r is taken to be the perihelion distance of the planetoid's orbit. 

An alternat«_ form of the temperature distribution, T(r) « r“^, has 
been suggested by Lewis (1974). In this case, mg “ r . 
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Once the critical mass is reached, a nucleus will accrete some 
gas along with the dust. As the mass increases still further, a 
larger fraction of the gas present near the nucleus wil'' be captured, 
so that in the limit of a very large mass the net density of 
captured material will be p = Kp^, which corresponds to the capture 
of all gas near the nucleus. For intermediate masses, the "effective 
density" of accreted matter is taken to be 

Pe “ KPd (K-l)]'‘, m > n^. 

which obeys the conditions p^ = p^ when m = m^ and p^ -*■ Kp^j when 
m -► » . This function is arbitrary and was selected by Dole primarily 
for its simplicity and its correct behavior in the limits. A 
functional form which is more physically exact would require 
knowledge of the structure of the accreting planet and a detailed 
dynamical analysis of the solar nebula, both of which are beyond 
the scope of this treatment. However, the expression is probably 
at least qualitatively correct and implies that the greater 
the mass of the pl^’-et. the greater the gas/dust ratio of the accreted 
mass . 

8. The nuclei are injected sequentially, with the newest nucleus 
growing to completion before the next is injected. The growth of 
the nucleus is calculated iteratively in the program, and "completion" 
is defined as a fractional mass Increase on a given iteration < 10”^ • 
Ideally, one would like to have all of the nuclei growing simultaneous! 
since, in the present form of the calculation, the final appearance 
of a planetary system is weakly dependent on the order in which the 
accretion nuclei are injected into the nebula. However, while this 
may change slightly the details of a given planetary system, the 
overall morphology of a set of planetary systems derived from 
similar initial conditions remains unchanged. 
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When the radius of capture of a growing planet Intrudes on 
that of an already -formed planet, the two coalesce into a new 
body which continues to grow until completion. The new semlmajor 
axis of the orbit of the component planet is taken to be 
a^ “ [m^^ + m 2 ] / [(mj^/a^) + (m2/a2)], 
where a^ and ag are the semimajor axes of the two coalescing bodies, 
and m^ and m 2 are their masses. The value a^ is the maximum 
allowed from the conservation of energy. The new eccentricity e^ 
is calculated from a^ and the conservation of angular momentum. 
Clearly, with no information about the position angles of the 
precollision semimajor axes of he two orbits, the three-body 
problem admits no unique solution and so (within the confines of 
the conservation laws) the choice of a^ and e^ is somewhat 
arbitrary. The form given above, however, is both physically 
realizable and convenient. 

9. A nucleus which is injected into a region which has 

already been swept free of dust by existing planets is a ’'dud*' 

and cannot grow, since a nucleus cannot initially accumul..te gas. 

Thus, the program ends when all dust between O .5 and 50 AU has 

been swept up. A typical run of the program will entail th^ 

injection of 100 - 300 nuclei, most of which are duds. The simulations 

in this paper were run on the IBM 370/168 at Cornell University. 

The running tim-- necessary for the formation of a single planetary 

system was on the order of 3 seconds, which (conveniently) is a 
15 

factor '^10 faster than the process being simulated. The cost 
was roughly $1 per solar system, or ten cents per planet. 
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III. THE ACCRETION PROCESS lUlNALi rAua 

« 

Most models of the solar nebula employ a self-gravltatinE 
disk or exocone < 1 AU thick and many AU in radius, with a total 
mass between 0.1 and 1.0 M^, 

exclusive of the mass of the Sun itself. Goldreich and Ward (1973) 

12 

have hypothesized a disk some 10 cm thick. As the nebula cools, 
the vapor pressures of some of the constituents fall below their 
partial pressures, and the condensation of small particles ensues. 

These particles then fall towards the central plane of the disk, 
accumulating matter as they fall from viscous drag and collisions in 
the medium. For 1 AU from the center of the disk, this occurs 
on a time scale 10 years and leads to particles with masses 
•v lOOg. This mass is an upper limit, however, being strongly 
dependent on the number of nucleation sites (i.e. the number of 
particles descending upon the central plane). ACRETE, .n injecting 
the nuclei sequentially, assumes a number of sites 100, as 
stated before. Hills (1973) suggests that there were 100 major 
accretion sites before mutual colllsJons led to fragmentation 

3 

into roughly 10 nuclei. The precise number, however, is of only 
marginal impoi-tance, for the resultant disk of particles in the 
central plane is gravitationally unst.ible and will clump together 
to form fewer pre-planetary accretion nuclei. This clumping leads 

to the formation of planetesimals with radii r 5km and masses 
1 fi 

m 'V' 10 g on a time scale of only a few thousand years. These 
planetesimals are largely in prograde orbits of near-zero inclination, 
and account plausibly for the coplanar nature of the solar system. 

In another model, Cameron (1973) states that turbulence in 
the solar nebula can cause grains to aggregate into bodies of a few 
tens of cm in radius which can then grow to lunar-sized planetesimals 
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as they descend to the central plane. As before, the process 
takes only a few thousand years. 

Once the preplanetary accretion nuclei have settled into 
the central plane of the nebula, their masses are much greater 
than the masses of the ambient dust particles, so that further 
growth will be dominated by the gravltatl'-nal capture mechanism. 
Weldenschllllng (197^) has performed a straightforward analysis 
of this process and concludes that, from accretion nuclei no 
larger than 10“^ of a terrestrial planetary mass, the solar system 

Q 

could be formed in about 10 years. 

The mass of the Injection nuclei which Dole used in ACRETE 

was m^ " 10“^^ Mq 'v 10^®gj coincidentally the same size as the 

planetesimals of Goldreich and Ward, but a good deal smaller than 

Cameron’s. We find that varying, the seed mass m^ by many orders of 

magnitude has absolutely no effect on the final results, since the 

amount of matter that the particle accretes from the nebula on the 

first Iteration is in most cases vastly greater than its initial 

mass. There., re, “ 10“^^ Mq was used in all subsequent runs. 

We can postulate a model similar to Goldreich and Ward's in which 

1 R 

numerous bodies of mass 10 g are created by local gravitational 
instabilities and which subsequently grow via accretion processes like 
those built into the computer program. Since the program is insensitive 
to the initial mass of the accretion nuclei, the number of nucleation 
sites in the solar nebula becomes unimportant. N’lmerous masses of 10^®g 
in Keplerlan orbits would eventually ccalesce into a s .aller number of 
more massive nuclei on which the accretion process would continue as 
before . 

Having established some theoretical basis for the specific model 



on which ‘program ACRETE is constructed, we proceed to alter the 
Individual parameters one at a tlme| to approach physically 
more realistic models than the one employed by Dole. 

IV. K: THE GAS TO DUST MASS RATIO 

In Figure 3 are displayed three model planetary systems 

generated with the canonical ACRETE program, but with different 

values of the ratio of gas to dust. As expected, decreasing the 

amount of gas In the cloud has no effect on those planets which 

never reach their critical mass m„ and hence never accumulate 

c 

any gas in the first place. This is apparent 
especially for the cases K = 30 and K = 10 in which the nuclei 
were injected into the same orbits in the two cases. (The orbits 
are determined by a random number generator: this is fed a seed 
number which causes the generation of a random series.). 

For K = 30 and K = 10, the five planets which did not 
accumulate gas (filled circles) underwent no change in mass, 
while the gas giants (open circles) are considerably smaller in the 
latter run. An extension of this result can be seen qualitatively 
in the K = 100 run, which has two very large gas giants. 

What is a reasonable value for K? Taking typical values for HI 
regions (Harwit, 1973), we find that the number density of grains 
is 'V' 10“^®cm“3, and that their radii are 'V' 3 x 10"^ cm. The mass 

-22 "3 

density of gas in an HI region is 2 x 10 g cm”J, so that if 
we assume unit mass density for each grain, we obtain K = 20. The 
value for a nebula of solar composition is 100, depending on the 
degree of condensation, but since K in the program is taken to be 
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the mass ratio of hydrogen and helium to all other substances (rather 
than the volatile: refractory ratio), this can be treated as an ^ 

upper limit. Thus, Dole's value of K * 50 is certainly an 
acceptable one, since even values as low as K = 10 in ACRETE 
yield plausible planetary systems. 

A more comprehensive model would have included the variation 
of K with r. Because the incidence of condensation should increase 
with declining temperature, K should decrease with heliocentric 

i: 

distance. However, we believe that a slowly varying K, or a 
bimodal distribution of K in which the values differ by a factor 
of no more than about 4, will not alter our results profoundly. 

We see from the figures that the only perceptible result of a 
variation of K by a factor of 5 is a change in the prevalence of 
gas giants. A distance-dependent K of the sort described would 
probably have as its principal consequence, a small inward 
displacement of the region of the Jovian planets. 
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V. THE CENTRAL DENSITY AND THE PARAMETER A 

Since An *»xponentlally decreasing density function leads to 
a total nebular mass which Is directly proportional to the density 

g 

at r » 0, changing the parameter A in the expression p^j « A exp(-ar ) 
is equivalent to scaling the mass of the cloud. More fundamental 
changes in the functional form itself will be discussed in a later 
section. This particular form was used by Dole because It has the 
mathematically desirable properties of being monotonically decreasing 
with r, and being integrable over a spherical or cylindrical volume; 
and because Its use in program ACRETE leads to the formation of planet 
ary systems resembling the solar system. To this latter end. Dole 
used the values o = 5 and 0 = 1/3, which for the moment we adopt. He 
employed the value 0.0015 Mq/AU^ for A, or roughly 10~^g cm"3, 
which leads to a total nebular mass of approximately 0.06 Mg when 
the opening angle of the exocone is taken to be tt/2 so that the 
"cone" is actually a sphere (see Appendix). (This mass, as the low 
central density indicates, is exclusive of the mass of the central 
star) . 

A mass of 0.06 M^ is somewhat low compared to that of most models 
Urey (197^), for example, derives a mass of 0.6 M^, although he refers 
to two other models which call for nebular masses of 0.2 and 0.05 M®. 
ACRETE does not "know" that some of 

the dust in the solar nebula is left unaccreted. In the actual 
formation process, it is possible that accretion onto planetary bodies 
from the solar nebula is in competition with a T Taurl solar wind 
tending to sweep away material. Workers concerned with the early 


stages of formation of the solar nebula often quote values of 
mass (exclusive of the mass of the central star) of about IMq 
( e.g. Cameron, 1976). But we are concerned with the values of 
solar nebula mass after the generation of accretion nuclei; 
in the interim a substantial loss of nebular material may have 
occurred, associated with the T Taurl stage of the central star. 

If we find that only a small rp-'ge of nebular density is consistent 
with familiar solar systems, it follows that such systems are 
correspondingly uncommon. The relative timing of the generation 
of accretion nuclei and the T Taurl stage of the central star 
is an important and as yet unresolved factor in understanding 
the origin of planetary systems. 

The results of decreasing the total mass of the cloud by one- 
third (A = 0.001 M„/Au 3) and by two-thirds (A = 0.0005 MVAU^) are shown 
in Fig. In the latter case, the density at every point in the 
nebula is so low that only one planet is able to accrete enough 
dust to exceed its critical mass and begin to accumulate gas. Only 
two planets have masses greater than one earth mass. For the case 
A = 0.001 M /AU the accretion process is somewhat more successful, 
although the resulting gas giants are small compared to those in the 
solar system. 

For values of A > 0.0015 M^^/AU^, Dole has presented some 
results (Fig. 5)^ Even a doubling of the total mass of the cloud 
leads to a near-catastrophic accumulation of gas by the large planets. 
Hydrogen thermonuclear reactions occur in the core of a star of mass 
^ 0.07 Mq, although deuterium burning will have set in long before 
that. Hence the planetary system generated by A = 0.003 M /AU^ in Fig. 5 
would quite likely be a borderline case of a double star system, while 
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those generated for A = 0.006 and A = 0.015 M®/Au3 would definitely be sc 
This is not necessarily a drawback to the model: the statistical 
studies of Abt and Levy (1976) suggest that virtually all stars are 
components of multiple systems, two-thirds of which Include 
stellar coirpanlons and one-third planetary. Hence the tendency for 
our accretion model to give rise to stellar or barely substellar 
companions to the central star mimics a similar tendency o nature, 
and the probability of the existence of numerous extrasolar planetary 
systems is correspondingly high. 

Thus .'or the numerical values and functional forms chosen, 

planetary systems of roughly familiar aspect are produced for 

nebular masses (exclusive of the central star) between about 0.02 

and about 0.2 M . f.ystems with smaller nebular masses than this 
0 

will tei'd to be comprised exclusively of terrestrial planets — 
and, eventually, of asteroids only. Systems with larger nebular 
masses w 111 evolve with the largest secondary components undergoing 
thernonucleai’ react ‘.ons, and therefore will become double or multiple 
star systems. In this case there will also be terrestrial and Jovian 
planets produced, some of which will be in orbits gravitationally 
unstable according to the restricted three body problem. But others 
will be in one of t> e three categories of reasonably stable orbits: 
around the center of mass of the systera if the two stellar components 
have a small reparation; around one or the other of the individual 
stars if ‘..e two components have a large separation; and in a figure-8 
trajecr.^ry around both components, though this is unstable in the long term. 

It is necessary to note, as Dole points out, that the generation 
of multiple : tar systems pushes program ACRETE somewhat beyond the 
limics its intended application (which for Dole was the simulation 



of planetary systems similar to our own). When planet formation 
gives way to star formation, ACRETE breaks down in the sense that 
the total mass of the companions can exceed the intended mass of 
the original nebula. This effect can be seen in Pig, 5d, in which 
the mass of the companions comes to 0.61 as compared with a 
nebular mass of 0.58 (derived from a central density of 0.015 

■3 

M^/AU and Dole’s radial density distribution). We will refer to 
such a breakdown of ACRETE as a pathological multiple star system. 

It arises from the breakdown of the approximation p(r-bp) = p(r+b^) 

(see section VII) when an accreting body becomes very massive. 

VI. THE ORBITAL ECCENTRICITY OP THE DUST PARTICLES: e 

In ACRETE, the accretion nuclei are assumed to capture all of 
those dust particles whose orbits cross their own. If a nucleus 
is injected with a high orbjtal eccentricity, it will, of course, 
cross the orbits of more dust particles, hence accumulate more of 
them and end up correspondingly more massive. Similarly, if the 
orbital eccentricities of the dust particles are high, then a given 
dust particle is more likely to cross the path of some nucleus. Thus, 
more eccentric particle orbits should give rise to more massive planets 
and, as a corollary, fewer planets in a given planetary system 
(since fewer nuclei are required to sweep up all of the dust). The 
results of varying e, the eccentricity of the dust particles in 
the solar nebula, between e = 0.1 and e = 0.5 are shown in Figures 6 
and 7. The results are as expected: a typical e = 0.5 run yielded 
six planets, three of which are quite large, while an e = 0.1 run 
yielded fifteen relatively small bodies. Curiously, for values of 

e ~ 0.3, the effect is nearly linear (Fig. 8; within the limits 
of uncertainty caused by the random injection of the nuclei (e = 0.1 



planetary systems will Generally have li|, 15, or 16 planets, etc.). 

The function, of course, must level off to N » 1, since for the 
limiting case e--*- 1.0 the dust particles will have near-parabolic 
orbits, all of which will cross the orbit of and hence be accreted 
onto the first nucleus Injected, leading to a double star system 
for all reasonable values of the nebular mass. 

We have argued in Section II, however, that frequent collisions 
in the early nebula would lead to a circularization of the orbits 
of accretion nuclei. This should apply to the dust particles as 
well, so that it is interesting that the program works "best "(in the sense 
of generating planetary systems similar to the solar system) when 
e = 0.25, which is a rather high value. Values of e more 
in accord with what we would expect in the nebula (say e s 0.1) 
lead to an inefficient accretion process when inserted into ACRETE 
(Pigs. 7b and 7c). 

We have attempted to counteract this effect by trying lower 
values of £ and increasing the mass of the cloud (by increasing the 
central density) to compensate. The results are shown in Figure 9* 
in which all three systems were generated with the same random number 
sequence. Figure 9a shows a run of ACRETE with Dole's parameters: 

A = 0.0015 Mq/AU^ and e = 0.25. In Figure 9b, A = 0.003 Mg,/Au 3 and e= 0.1, 
so that the mass of the cloud is now * 0.12 sin ©^ax condition 

was shown earlier to give rise to a system of barely substellar companions 
when £ = 0.25, as shown in Figure 5b. Now the accretion process has 
not run away as dramatically; the largest gas giant is only seven times 
the mass of Jupiter. 

The effect of lowering e still further by another factor of 
five, to e = 0.02, and Increasing A again by only 25 % (to A = 0.00375 
^©/AU^) is illustrated in Figure 9c. The mass of the largest companion ha 
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increased by 50 Jf, to 0.01 M^. Clearly, we have pushed A almost 
to the limit; decreasing e to zero and increasing A much further 
will lead to a binary star system. 

With nebular masses of the order of 0.1 the eccentricity 
of dust orbits of roughly 0.15 seems to produce recognizable 
planetary systems. Both values seem to be in reasonable conformity 
with our expectations for the solar nebula. We also note that 
the mean eccentricity of the asteroids is e = 0.15. However, 
because the accretion process as simulated in ACRETE produces 

■3 

familiar solar systems when e = 0.25 for A = 0.0015 M^/AU , 
subsequent computer runs will for convenience use that value, 
as we vary other parameters. 

VII. MODIFICATIONS TO THE DENSITY DISTRIBUTION 

The mass density as a function of heliocentric distance is 
a critical attribute of any model of the solar nebula. The functional 
form of the density which Dole used in ACRETE is = A expC-ar^*^^) , 
with A and a as free parameters. More generally, we may use the form 
(which we will call form A) = A exp(-ar^), and treat 3 as a free 
parameter as well. The mass which a nucleus at radial distance r will 
accumulate is roughly proportional to r^p. For form A, with 3 = 1/3, 
this reaches a maximum when the r derivative of r“^p^ vanishes; 

i.e., when r^^^ = (9/a) . Dole used the value a = 5, which leads 
to the. largest planets near r = 5*8 AU, i.e. 

not very far from Jupiter’s orbit. Thus, we can move the position 
of the largest planets by altering the value of a. Large values of a 
will make the exponential drop off faster, so that distant planets 
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3 

become smaller as the maximum of r p moves Inward. 

The assumption of a solar nebula with a density maximum at 
10 AU Is supported by the frequency histogram of separations 
of double star systems (Kulper, 1951) which Is also peaked near 
10 AU. This Is In reasonable accord with the Implications of 
prograun ACRE3E, that In many cases the formation of double stars 
Is due to the condensation of a particularly massive Jovian planet 
from a solar nebula. If we are to preserve the total mass of the 
nebula as a Is Increased we must Increase the value of A. If we take 
an exocone of angle 0 = ir/2 (l.e. a sphere), then for a density 
distribution with form A and $ = 1/3, the total mass of the cloud is 
(see Appendix) = ^BSS^OTiKA/a^MQ. For K = 50, A = 0.0015 M^/AU^, 
and a = 5, this becomes = 0.06 If we keep K = 50 

and wish to preserve « 0.06 M^, then the relation between A and a 
Is A = 0^/(1. 3 X 10 "^) M^/AU^. The results of varying a (and A with a) 
2U7e shown in Figure 10. For small values of a, more planets would 
be formed at r >> 50 AU if the program were allowed to inject 
accretion nuclei out that far. For a = 1, for example, Jovian planets 
would be formed near 700 AU. 

It Is apparent that, although In principle an exponentially 

decreasing density distribution is reasonable, the particular form 
1/3 

exp(-ar ) is quite arbitrary and only serves well for values of a 
not very different from 5. This form was used in the first place 
because of its pleasing tendency to produce familiar end results. 

In fact, however, it falls off much more rapidly than most other 
theoretical models, decreasing to 1 percent of the central density 
at only O .78 AU (Just outside the orbit of Venus). This indicates that, 
for such a model, planetary formation takes place at a very late stage 
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In the condensation of tlio nebula. 

Heppenhelraer (197*0 has manaccd to find a meetinc point 

between Cameron and Pine's models and a density distribution of form A, 

by matching the pressure Implied by such a distribution to the pressure 

obtained by Cameron and Pine (their Figure 2). If we assume that the 

nebular material obeys the ideal gas law Larson's (1969) adiabat 

gives p « Hence Heppenheimer fits the relation P = Pq exp[-(5/3)ar^ 

to the pressure curve of Cameron and Pine and finds a = and 6 = 0.22. 

The mass of the nebula then becomes 

M = 1.18 X 10^ A sin 0--, M 
c inax e 

In the model of Cameron and Pine, the nebula has rcnghly the exocone 
geometry, with a semithickness of approximately 1 AU at about 50 AU 
from the center. Hence = 0.02, leading to « 236A for A 

3 

given in M /AU . If, to retain consistency with other models, we 

demand = 0.1 M_, we find A = *1.2 x 10“^ M /AU^. 

Running ACRETE with the parameters a= *<.*], 3 = 0.22, and A = i|.2 

X 10”^ leads to a pathological multiple star system. Clearly, this 

3 

same result will occur if A is increased to 0.0015 M@/AU , in which 
case all of the parameters would be identical to those used by Dole 
except for 6, which is 0.22 instead of 0.33- The companion star, then, 
forms near the edge of the nebula; Table 1 indicates that it is the 
more gradual decrease in the former case that is causing the difficulty. 



TABLE 1 


DEIENDENCE OP DENSITY ON HELIOCENTRIC DISTANCE FOR TWO MODELS OP FORM A 


r exp(-5r^^^) 

0.1 AU 9.82 X 10“2 

1 6.73 X 10"3 

2 1.84 X 10“3 

10 2.10 X 10"5 

20 1.28 X 10"^ 

50 1.00 X 10"® 


exp(-4.4rQ«^^) 
7.06 X 10“2 
1.23 X 10“2 
5.9‘J X 10“^ 
6.74 X 10“^ 

2.02 X 10"^ 

3.03 X 10"5 


It is nonetheless striking that a fairly small change in the 
qualitative behavior of p should lead to such violent changes in the 
results, especially in light of the fact that the form exp(-4. , 
which is even more similar to Heppenheimer's result, gives rise to 
a comfortably familiar planetary system (see Figure lOd). 

A more complicated aspect of the entire problem — and one which 
was not taken into account at all in the original version of ACRETE — 
is the variation of nebular density with z, the distance perpendicular 
to the plane of symmetry. Since in general the protoplanetary orbits 
will have some nonzero inclination to the central plane, the mass 
accreted on each orbit will then be determined by the nebular surface 
density at appropriate heliocentric distance. ACRETE, however, assumes 
orbits with zero inclination so that the only accreted material is 
that within the toroid (shown in cross section in Figure 11) defined 
by the orbital eccentricities of the accretion nuclei and dust 
particles and the gravitational cross section of the nucleus. The 
volume of the toroid is approximately V » 2Trr(bg^ + hp)(x^ + Xp), where 
Xg and Xp are the gravitational capture distances of aphelion and 
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perihelion (see section II), and b and b include the effects of 
the dust’s orbital eccentricity. If the density does not vary 
out of the plane, and if p(r - bp) (r + b^), the mass of dust 
within the toroid is - 2Trr(b^ + bp)(x^ + Xp)pjj(r), a formulation 
employed by Dole. 

From studies of a self-gravitating set of mass points, 
originally due to Ledoux and to Chandrasekhar, Urey (197’0 argues 
that the vertical variation takes the form 

p(r,z) = p(r) sech^ [z/H(r)], 

where H(r) is some radlally-dependent characteristic vertical 
scale height. In this case, the mass of the dust within the toroid 
becomes f ] 


M- * 2'Trr(b. 


= 27Tr(ba + b„)H(r) 


I'p) I sech2 dz + l^p sech^ 


dz 




tanh [x^/H(r)] + tanh [Xp/H(r)] 


Pd(r,0) 


Pd(r,0) 


VH(r) ->■ 0. Urey 


a • -p' 

\ 

which approaches Dole's form in the limit Xj, 
finds (his Table III) that H(r) = 0.0026?r AU; for r measured in AU. 
For a planet of mass m with a circular orbit, we have x^^p/H - 
374(m/M^)^^\ so that tanh (x./H) will differ appreciably from its 
argument (and therefore depart from Dole’s limiting case) whenever 
m K Since this is < 10”^ the mass of the Moon, we can 

conclude that, for every case of Interest, Urey’s vertical density 
distribution will lead to results substantially different from Dole’s 
vertically-unlform model. 
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A sech 2 vertical decrease In density Is faster than exponential. 
Urey's nebula, therefore, Is much thinner than Dole's. A characteristic 
scale height of 0.00267r leads to a semithickness of roughly 0.1 AU 
at r “ 50 AU, or about an order of magnitude thinner than the nebular 
model of Cameron and Pine. Hills (1973) also concludes that the 
scale height for gas and dust in the solar nebula is very small: roughly 
0.1 r AU for Hg gas and lO”^ r AU for dust. For models which are :,o 
concentrated into the central plane, however, radial density ■ rl- 
butlons as steep as Dole's (see Table I) lead to very low nebi. 
masses. Consider an exocone of radius R and opening angle 6 

max 

In which the only density variation Is radial, p(r). Then the mass 
of the cloud is 

”o,Dole ‘ Vi 

^0 

Urey's nebula has a vertical density distribution as well, and has 
cylindrical geometry which we will characterize by some thickness h. 

Then the mass becomes 

fR p 

^c,Urey Lh sech^(z/yr) drdz 

where yr - 0.00267r = H(r). If we let h >> yr (true for h > 0.25 AU) 

we can, to good approximation, extend the limits on the z integral to 

+ 06 and Integrate to get 

Me, Urey “ yr2p(r)dr 

^0 

so that ^c,Urey/^c,Dole “ Y/sln 6n,ax’ Cameron and Pine's nrbula, 

®max “ 0.02, so that if we were to apply this opening angle to Dole's 
exocone, we could get ^e ^urey/^c,Dole a complete sphere, 

the ratio becomes 2 x lO"^. Thus, for Identical radial density 
distributions, Urey's model leads to nebular masses much smaller 
than Dole's. 
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In principle the nasaes of the two models could be recon- 
ciled by either Increasing the central density in Urey’s mod»-l 
or inserting a radial density function that is less steep than 
Dole’s. We have simulated a sech (z/yr) dependence in the vertical 
direction by modifying the mass contained in a toroidal volume in the 
fashion derived earlier. Using a radial density function of form A 

with ot = 5* 3 = 1/3 (Dole’s values) and a central density increased 

•a 

by an order of magnitude from Dole’s value (to 0.015 M^/AU ) , ACRETE 
consistently generates pathological multiple star systems. Nonpatho-? 
logical systems are generated wich any degree of regularity only 

O 

when the central density becomes Po li 0.006 M^/AU . This implies 

”c,Dole " “-23 sin M, 


^CjUrey ~ 10 ^ 


M 


which is unrealistically low for reasonable values of 6. The 
failure of ACRETE to accommodate a vertical density distribution 
must unfortunately be interpreted as being due to the exclusion of 
some pertinent physics from the computer model. It seems clear 
that a vertical distribution will exist in a rotating preplanetary 
nebula and that th« accretion nuclei, like the planets, will travel 
in orbits with nonzero inclination. The Importance of the surface 
density (as opposed to the vertical and radial volume densities) 
arises from these conditions, but is overlooked in the computer 
program in the assumption of pe fectly coplanar orbit:',. 


* ;■ 
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VIII. THE ARBITRARY NATURE OF p(r) 

A rigorous derivation of the radial density distribution 
requires a detailed knowledge of the equation of state at all 
points in the cloud, a treatment which is well beyond the 
scope of this paper and, apparently, many others. Numerous 
assumptions, including consideration of the ambient magnetic 
field and the solar wind enter into the problem, and the 
final results mu.t be strongly model-dependent. 

Any formulation of the density function in the cloud 
must at present be, to a certain extent, arbitrary, so that 
it is perhaps the safest course to choose one which is charac- 
terized solely by physically reasonable qualitative attributes. 
This is essentially what Urey, Dole and Cameron and Pine all 
did, and what we shall proceed to do. 

One of the most obvious forms to try, because of its 
simplicity and wide applicability, is a simple exponential. 

This is just another manifestation of form A, with B = 1. 

g 

Since form A contains r , we see that, because 1 > 1/3, P will 
fall off more quickly than Dole's form. Various radial scale 
lengths and central densities were tried, with the most 
success in generating planetary systems obtained with scale 
lengths of 0.4 + 0.1 AU. The results of the density 
1 unctions 

p(r,z) = 0.004 e~^^ sech^ (z/yr) M /AU^ 

© 
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p(r,z) = 0.002 sech^ (z/yr) M^/AU^ 

(for r in AU) are shown in Figures 12b and 12c. The masses of the nebulae 
in these cases are 10 and are kept deliberately low because 

of a propensity for this form of the density distribution to generate 
pathological multiple star systems. Removing the vertical density 
dependence does not help; the resultant increase in mass concentrated 
into the inner region of the nebula only exaggerates the tendency 
towards pathological results. 

Another obvious form to try is a power law, p(r) « r“” where 
n > 0. This form has the disadvantage .of diverging at zero, although 
this is clearly not a problem physically, since we are only interested 
in the nebula at r i 0.1 AU. For the form to be integrable, we further 
require that n > 2 in a strictly cylindrical nebula, and n > 3 for a 
spherical one (although this is not a rigid restriction since the nebula 
has a finite diameter) . We can deal in another way with the divergence 
at small r with the following ad hoc argument: most models of the 
formation of the solar system suggest that the young sun was in a T Tauri 
stage during the epoch of planetary formation. The T Tauri solar wind 
would considerably deplete the interior portion of the nebula of both 
refractory and gaseous material. Evidence for the size of a depleted 
region is suggested by recent observations of the T Tauri star RU 
Lupi by Gahm, et al. (1975) > who found concentrations of dust, 
presumably driven out by the stellar wind, at distances of a few 
tenths of an AU from the star. We therefore choose to modify our 
power law distribution so that it reflects some flattening of 
the mass density function in regions close to the star. This 
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we refer to as form b: p(r) ® Pj^(r” + C)”^ vfhere C Is some 

dimensionless constant and r Is measured In AU. p(r) approaches 
the simple power law r’*’^ when r >> For a given n, we can 

solve for p^^ and C by demanding a particular central density p^ « p-^/C 
and a particular nebular mass, which for Dole’s exocone Is 
»o,Dole ■ "’"'Pi I" Pin M, 

If Pj^/C Is the central value of P^(r) and R Is the radius. In the 

symmetry plane, of the nebula. Prom the arguments given earlier, 

2 

when one includes the sech z vertical distribution the mass becomes 
approximately M^^urey ' ”c,Dole 7/sln 0 .Tjax, where y = 0.0026?. M^, 

Is evaluated for various values of n In the Appendix. 

Note that for any density distribution of form B, the function 
p3p reaches a maximum at r^j^ = (3C/n-3)^^"^. As stated earlier, rj^^ 
represents the distance at which the la*rgest planets In the system 
will tend to form. 

If we choose to simulate Urey's r“^ density function with 

the T Tauri modification, we must decide on R, since, with an Infinite 
2 3 1 

upper limit, r (r +C)“ is not integrable. Taking R = 70 AU and, as 
before, K = 50, the mass becomes 

“c.Dole * 67i.p^ln{l+70 C'1/3) sin 6 

Letting = tt/ 2, the conditions = 0.06 and p^ = p^/C = 

0.0015 M /AU^ yield the approximate solution 
© 

P.(r) = ^ M /AU^ 

^ + 0.032 ® 

The distribution remains fairly flat out to r 0.032^^^ = 0.3 AU. 

A planetary system generated by ACRETE with this density distribution 
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ls shown In Figure 13a. Two of the planets are very large, but are 

still sub-stellar (the largest Is ten times the mass of Juriter) . 

Note, however, that for n ® 3* ^as no maximum at a finite r : 

m* 

In models of this type, the largest planets will always be formed 
at the outer edge of the nebula. 

2 

If we Insert the sech z vertical distribution, a tendency 
towards pathological multiple star systems develops. As before, 
we can combat this only by lowering the total mass of the cloud by 
a factor that Inhibits the formation of gas giants. The resultant 
planetary systems, generated when urey resemble the 

one In Figure 10a. 

The choice of n is arbitrary and is open to considerable 
experimentation. When n = 6, for example, the mass Integral yields 
the result 

.2 


”c,Dole “ (2irV3)(Kp;^C-l/2) gin 6 


max 

In this case, the Integral converges as Taking Smax = tr/2, 

M.,,Doie = 0.06 and Po = p^/C = 0.0015 M^AU^ as in the n = 3 case, 

the density becomes 

M /AU^ 

^ ~ r& + 0.01^1 

Now r^p reaches a maximum at r^^^ = (0. 052/3 0.5 AU, out to which 
distance p is flat. The result of this distribution is shown in 
Figure 13b. Gas giants can only form close to the sun, followed 
by terrestrial planets and asteroids, moving outward. That such 
a planetary system can form at all is highly questionable; the 
inner terrestrial planets might not have stable orbits, and both 


r 
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the T Taurl wind and Jeans escape would make the accretion of 
larco amounts of c^s by a planet so close to the sun unlikely. 

Introducing the vertical density function causes the same 
problems as before. A low cloud mass permits the formation 
of terrestrial planets close to the sun, but otherwise gas giants 
turn into stellar companions. 

Finally, we note that Larson's adlab.at suggests a pov.'er law. 
If T(r) « and we use the simple temperature law mentioned 

earlier, T(r) « r"'*, we obtain p(r) « If T(r) « r“^, 

p(r) « Such extremely shallow distributions, however, 

generate only pathological systems. 


IX. PLANETARY DISTANCES 

The geometric spacing of the planets in the solar system Is 
one of its most striking properties, represented by a number of 
formal schemes the most famous of which is the so-called Titius-Bode 
"law" — in which the semi-major axes of planetary orbits in AU are 
written r = 0.4 + 0.3 x 2^. The value -«must be assigned to n 
in order to explain Mercury; thereafter, integer values are adopted 
beginning with zero. It is then necessary to identify Ceres as a 
• lanet. Even so the values for n - 7 (Neptune) and n » 8 (Pluto) 
are in unsatisfactory agreement. Thus Bode's law can be described as 
a fit to eight numbers by an equation with five or six free parameters 
or arbitrary indexing conventions — not a very impressive "law". 

Derraott (I 968 ), however, has proposed a simple, quasi-geometric 
form which describes adequately the spacing of satellites around some 
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of the major planets and meets wluh moderate success when applied 

to the solar system as a whole. If Is taken to be a constant 

of proportionality then the periods of the planets can be expressed 

approximately as ® where J is a small integer (J = 6 

for the solar system) and n is a given planet's "orbital" Integer, 

generally about the same as its serial position outward from the sun. 

Differences between n and the serial position arise because Dermott 

allows that two planets can share the same value of n and furthermore 

that all values of n in a sequence need not be used. For the solar 

system, both Earth and Venus are in the n = 2 orbital, and both 

Neptune and Pluto share the value n = 8. The advantage of this "law" 

over Bode's is that the relationship between P^ and n can be graphed 

as a straight line semilogaritlimically j in light of the amount of 

freedom in the choice P , J, and the n's, however, it is probably no 

c 

less arbitrary. 

A measure of the adequacy of the relation, used in part by Dermott, 

can be made by comparing the n's to the m's in the equation P_ = 

here, the Pj^^ values represent the actual periods of the planets, and m 

takes on any values (not necessarily integer) to ensure that that is the 

case. Taking Dermott 's quantity An = m-n, we define the quantity 

a = [Np^ I (An)^]^*^^: the rms derivation from the law per planet when 

Np is the number of planets in the system. 

In order to compare our simulated planetary systems against 

3 

Dermott 's law, we reformulate the latter as a = C J , using Kepler's 

n a 

third law to utilize the orbital semimajor axes a^ rather than the period 
(all appropriate constants are now absorbed into C^^). Values of j, C^, 

Np and a for the solar system and for some of the systems generated 
in this paper are shown in Table 2. Note that allowing a half-integer 
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value of J (6.5) for the solar system results in a significant 
inq>rovement in o over the integral case; this serves to illustrate 
the somewhat arbitrary nature of the procedure. It must also be 
borne in mind that even for randomly distributed values of m and n, 
the rms value of A n is I/ 2/3 “ O. 289 . Since o is nearly half this 
value in even the best case^ it is apparent that our model planetary 
systems follow a Bode-type law about as well as the solar system. 

The agreement of Bode-type laws with our model solar systems 
even in such bizarre cases as, say, 12a or 13a is of some interest. 

It cannot be due to multiple resonances in the n-body problem as 
proposed by Molchanov ( 1968 ) because the appropriate physics is not 
contained in the computer simulations [See also other criticisms and 
Molchanov’s reply: Backus (1969), Henon (1969)> Molchanov (1969), 
Molchanov (1969)* Gingerich (1969), Dermott (1969)3. 

Instead, what seems clearly to be happening is a kind of collisional 
natural selection. The solar system begins with gas, accretion nuclei and 
dust grains, and a variety of orbital eccentricities and heliocentric 
distances. But because of the high sticking efficiency in nucleus- 
grain and nucleus-nucleus collisions, those accreting planets with 
Interacting orbits merge. In all cases the final configuration shows 
planets nicely separated one from another. Because larger quantities 
of mass are required to generate the jovian planets, they are required 
to sweep up larger volumes of dust and therefore have larger mutual 
separations than do the terrestrial planets. Because there were then 
more objects on more eccentric orbits in a time before the completion 
of this collisional natural selection, the rate of planetary collision 
very early in the history of the solar system may have been 


considerable — quite apart from the infall of matter in debris 
rlnES in the vicinity of forming planets. 

TABLE 2 

PITS OP REAL AND MODEL SOLAR SYSTEMS TO A MODIFIED DERMOTT RELATION 


System 

N 

-P . 

(AU) 

i 

0 

Comments 

Solar System 

10 

0.263 

6 

0.217 

Includes Ceres 

Solar System 

10 

0.230 

6.5 

0.169 

22 X Improvement 
in 0 

Figure 6(a) 

6 

0.185 

8 

0.195 

Form A 

Figure 12(a) 

11 

0.236 

4 

0.174 

p 

sech z vertical 
distribution 

Figure 13(a) 
1 

8 

0.235 

7 

0.134 

p(r) “ l/(r^+C) 


X. CONCLUSIONS 

The original results of Dole’s program are so provocative that 
it is natural to question whether the remarkable planetary systems 
generated by ACRETE are the result of careful tailoring of the 
assumed radial and vertical density distributions in the solar nebula 
and the accompanying free parameters; or whether they are properties 
of any reasonable set of assumptions about the solar nebula. We have 
confirmed that for a perhaps plausible radial density law (form A) 
add the arbitrary choices a = 5 and B = 0.33 recognizable planetary 
systems of solar system type are generated for nebular masses 
between 0.02 and about 0.2 Mg and dust grain orbital eccentricities 
not extremely high or extremely low. However, in Section VII we found 


























that an apparently small change from B ■ 0.33 to B ■ 0.22 leads 
to a striking change In the end product. It is, of course, possible 
that a correct reconstruction of the underlying physics of the solar 
nebula will yield values of B near 0.33 and density distributions 
like form A. But in the absence of such a justification we can only 
conclude either a) that ACRETE is missing some of the essential 
physics of solar system cosmogony, or b) that planetary systems 
of our type are only one example in a rich array of alternative 
varieties of planetzucy systems. Likewise, more fundamental changes 
in the nebular morphology (e.g. from an exponential to a power law 
density distribution function) generate planetary systems some of which, 
although they do not closely resemble our own, are not fundamentally 
objectionable (Figs. 12 and 13). 

Abt and Levy (1976) have found that the frequency of secondary 
masses for binaries with periods less than a century varies as the 
one-third power of the secondary mass. If this function can be 
extrapolated, it implies that about 20 percent of stars of solar mass 
have a largest companion of mass 'v 10 “^Mq and about 10 percent of stars c 
solar mass have a largest companion of mass IO^'^Mq. In at least a crude 
way this result is consistent with our findings: where the solar 
nebular mass is between about 1 and 10“^ binary stars form; 
while for smaller solar nebular masses during accretion, jovian 
planets form. (This is a model-dependent consistency, however, 
because Abt and Levy believe that the short period binaries are 
fission systems from a single protostar). 

Perhaps the most striking result of tnls exercise is that in 
all cases planetary systems are generated which satisfy a Titlus- 
Bode sort of law; and in all cases the number of planets generated 



- 3*1 - 


is between several and about 20 — that is 10. This property 
may be understood in a very qualitative way: an accreting planet 
can be expected to perturb the orbits of dust particles in the 
solar nebula up to a few AU distant, depending on the planetary 
mass, and to sweep up material in such a zone. The number of 
such zones of a few AU in width in a nebula 50 AU in radius is 
V 10; hence the number of planets. The remarkable result on 
the number of planets then is attributable to the size of the solar 
nebula which is given as an input parameter. However, there is at 
least a hint in our results (see Pig. 13a) that more massive and 
extensive solar nebulae lead to very large jovlan planets or very 
small stars; but that, even in such a case, the number of dust lanes 
swept up is still 10. The width of the lanes, i.e., the 
spacing of the planets, arises in part from the dependence of 
the accreting planets* gravitational capture lengths on their 
orbital radii, x « r (see Section II). 

The computer simulations described in this paper and in 
Dole’s take no explicit account of chemical fractionation, the T 
Taurl stage of early stellar evolution, nebular opacity, frozen- 
in magnetic fields and a number of other factors. The accretion 
process is imagined to be purely dynamical, and in that respect 
is similar to the work of Weidenschilllng (197^). Furthermore, 
for purposes of compucatlonal convenience even the dynamics is 
simplified; as, for example, when all dust grains are taken to 
have the same eccentricity. We have found that both major and 
minor changes in the model, of equal apparent plausibility as 
the initial conditions assumed by Dole, lead to dramatic changes 
in the resulting planetary systems. Of course, it must be borne 



- 35 - 


rtproducbility of the 

PAGE IS POOR 


in H!ind that with at least six free parameters to describe the 
model nebula, we have investigated only a small fraction of our 
’’parameter space." In fact, there exist a multitude of such 
spaces, each defined by distinct, plausible density distribu- 
tions, of which only two — forms A and B — were considered in 
this work. Considering these, we have covered a relatively 
minute number of cases indeed. 

Nevertheless, in all cases — even when pathological binary 
star systems are generated — planetary systems are formed. One 
Interesting result is that while terrestrial planets can be 
formed without Jovian planets (in very low mass solar nebulae) 
the converse never occurs. We continue to be impressed that so 
simple a dynamical model generates recognizable if not familiar 
planetary systems with 'v ten planets per system and a Bode’s law 
spacing for a wide variety of initial conditions. The results 
suggest that planetary systems are widely prevalent in the Milky 
Way Galaxy, but that substantial morphological differences between 
extrasolar planetary systems and our own can be expected. 
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APPENDIX: EVALUATION OP THE MASS OP THE NEBULA 


The total mass density throughout the cloud (both gas and 


dust) Is 


p(r) “ P^(r) “ Pd ^Pd ~ Kp 


since K 50 >> 1. Thus, for the geometry of an exocone of opening 
angle and radius R, the mass of the cloud is 


M = ^irK 
c 


Pd^ 

0 


r)r^dr sin 6 
' max 


For a cylindrical geometry with Urey's sech^z vertical distribution 

this becomes (See Section VII) 

fR o 

M = 4itKy p^(r)r‘^dr 


g 

If the density distribution takes form A, = ke~^^ , the mass in 


fn p 

M = 4 KA5 r^e dr 
A 


either case is 


where C ® sin Qj^ax ^ exocone, and C = Y = 0.0026? for Urey's 
cylinder. Typically, R 'V' 50 and a'V'5. If 6 = 1, the centroid 
of the distribution is at r^ = 1/5 « 50, so that, to excellent 
approximation, we can extend the upper limit on the the integral to 
“ . This approximation is still quite good for 6 <1, so that we 
can make the substitution u = r*^ and perfor... the semi-infinite integral 


to find 


r (3, 


If the distribution takes form B, = 0 ]^(r^ + C)"^, the integral 


'R r«2(jT» 

”b ' ^ 

. 0 


iiTK^^tan-^^ n = 6 

3 /TT /c 
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as R 


The approximation R • is good provided only 
1/2 

condition R/C » 1 holds, and, since R 'v* 50 
is generally true. 

If n » 3* does not converge as R . 
Integration is 

Mr “ ^ KCPi *in(l + RC"^/3) 


that the weaker 
.ind C 'X' lO’^, this 

The result of the 

n - 3 


Fortunately, the divergence Is logarithmically slow, so that the 
choice of R is not c.-ltical. Since ACRETE is usually run so as to 
inject nuclei out as far as r = 50 AU, we calculate Mg with R = "’0. 
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FIGURE CAPTIONS 



i. 


Pig. 1 


Pig. 2 
Pig. 3 
Pig. 4 
Pig. 5 

Pig. 6 

Fig. 7 
Fig. 8 

Fig. 9 


Planetary systems generated by Dole (1970) using ■ 

program ACRETE. Solid circles represent terrestrial 
planets, while open circles indicate Jovian planets - 

that have accreted gas as well as dust. The radius 
of each circle is scaled solely by the cube root of , 

the planet's mass^ given- in the figures in units of 
Earth masses. The positions of jovian planets are given 
by the centers of the circles. Those cases in which 
planets to be abutting or overlapping are, of course, - 

only artifacts of the schematic diagram. All planets 
produced are well separated, as the relative positions * 

of their centers indicate. The fourth system displayed , 

is our own; the others are generated by ACRETE. I 

The exocone. In the original model, the density of ' 

gas and dust depended only on the radial distance ! 

from the central star. 

Planetary systems generated by ACRETE for different 

values of the gas: dust mass ratio in the solar nebula. | 

Dole employed the value K = 50. - 

The effect of decreasing the central density c the ' 

nebula, which is equivalent to scaling the mass of ^ - 

the cloud. Densities are measured in units of Mg/AU^. • 

The effect of increasing the central density of the 
cloud from Dole's value of 0.0015 Mg/AU3. In the patho- 
logical system (d), the sum of the masses of the bodies j 

slightly exceeds the original mass of the nebula. h 

Increasing the eccentricity, e, of the orbits of the 
dust particles in the nebula results in a more efficient 
accretion process and a higher incidence of Jovian 
planets. Results which resemble the solar system are 
obtained when e = 0.25 as shown in (c). 

Decreasing e leads to numerous terrestrial bodies and ; 

relatively few gas giants. Small values of e may have s 

prevailed in the actual solar nebula. 

i ' 

Average number of planets N generated by ACRETE as a j' 

function of e, the orbital eccentricity of dust particles 
in the solar nebula. For e i 0.3, the effect is nearly If 
linear. 

I 

Compensation for the inefficiency of the accretion process i 
associated with small values of e by increasing the mass 
of the nebula, proportional to A. Although the results ;; 

illustrated in (b) compare favorably with the system = 

shown in Fig. 5(b), the central density A cannot be 
increased much beyond the value 0.00375 Mq/AU^, as shown 
in (c), without generating binary star systems. n 







Fig. 10 


Pig. 11 


Fig. 12 


Pig. 13 


Systems generated by ACREl. .y varying the steepness 
with which the density prof-ie falls off. a is the_ 
parameter in the expression p, = 0.0015 exp (-ar^'^). 

In cases (a), (b) and (c), more planets with orbital 
semimajor axes r > 50 would be formed if injected 
accretion nuclei were allowed out that far. 

Cross-section of the toroidal volume swept out by an 
accretion nucleus with orbital semimajor axis r. 

X and x are related to the gravitational cross-section 
of the nucleus, b^ and b are also related to this 
cross-section as well as to the orbital eccentricities 
of both the nucleus and the dust particles in the nebula. 

2 

The effect of including Urey's sech z vertical density 
profile. Figure .£a) was generated by the insertion of 
0.0015 exp (-5r^'^). Figures (b) and (c) were obtained 
by using a simple exponential radial density profile, 
with scale lengths of 1/2 AU and 1/3 AU, respectively. 

Planetary systems obtained with radial density distri- 
butions different from an exponential. ^ 

Figure (a) was derived for the form p. = p,(r^ + C 
Figure (b)for P^j = + C)“^, where p. and C are 

constants. Figure (c; shows the solar system. 
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